Spontaneous Symmetry Breaking in Directed Percolation with Many Colors: 
Differentiation of Species in the Gribov Process 
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A general field theoretic model of directed percolation with many colors that is equivalent to a 
population model (Gribov process) with many species near their extinction thresholds is presented. 
It is shown that the multicritical behavior is always described by the well known exponents of 
Reggeon field theory. In addition this universal model shows an instability that leads in general to 
a total asymmetry between each pair of species of a cooperative society. 
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It is well known that the Gribov process yJJ^ (Reggeon 
field theory |^|^, stochastic Schlogl model is a 
stochastic multiparticle process which describes the es- 
sential features of a vast number of growth phenomena of 
populations near their extinction threshold and without 
exploitation of the environment. The transition between 
survival and extinction of the population is a nonequi- 
librium continuous phase transition phenomenon and is 
characterized by universal scaling laws. The Gribov pro- 
cess belongs to the universality class of local growth pro- 
cesses with absorbing states such as directed per- 
colation [p|-pl|, the contact process p^p^, and certain 
cellular automata ||l^,|l^ , and is relevant to a vast range 
of models in physics, chemistry, biology, and sociology. 
I have been studying this transition for some time by 
renormalized field theoretic methods 0|. In this letter I 
will show that the Gribov process with many interact- 
ing species (directed percolation with many colors) also 
belongs to this universality class. In addition this mul- 
ticolored directed percolation exhibits a new interesting 
phenomenon: even if the dynamics of the different col- 
ors is modeled symmetrically, there exists an instability 
which leads to a differentiation of the species in their ac- 
tive state. This color symmetry breaking leads to many 
different stable forms of a cooperative society. 

The multiparticle processes under consideration are 
characterized by the following four principles: 

1. Self- reproduction ("birth") is error free, but may 
be different for each species; particles can be 
annihilated with rates depending on the species 
("death"). 

2. The particles interact with species-dependent cou- 
plings ("competition" and "cooperation"). 

3. Particles diffuse in a d-dimensional space ("mo- 
tion"). 

4. The states where some species are extinct are ab- 
sorbing. 



In the language of chemistry, the system corresponds to 
an autocatalytic reaction scheme of the form Xa <-> 2Xa, 
Xa — * 0, Xa + Xp kXa + IX whcrc the last reac- 
tion subsumes the interactions of species with different 
colors a, /3, and k^l — 0,1. A description in terms of 
particle densities typically arises from a coarse-graining 
procedure where a large number of microscopic degrees of 
freedom are averaged out. The influence of these is sim- 
ply modeled as Gaussian noise-terms in a Langevin equa- 
tion which have to respect the absorbing state condition. 
Stochastic reaction-diffusion equations for the particle 
densities n (r, t) — {ni (r, t) , 712 (r, t) , . . .) in accordance 
with the four principles given above are constructed as 

dtUa (r, t) = XaV'^Ua (r, t) + Ra {n (r, t)) (r, t) 

Ha{r,t) , (1) 

where the first term on the right hand side models the 
(diffusive) motion, and the are the reproduction rates 
of the particles with color a. These deterministic terms 
are constructed proportional to tiq in order to ensure the 
existence of an absorbing state for each species. Expand- 
ing the rates Ra in powers of n results in 



Ra {n (r, i)) = -A^ 



(2) 



up to subleading terms. The Gaussian noises Cq (r^ t) 
must also respect the absorbing state condition, whence 



(Ca (r, t) C/3 (r', t')) = \aga5a,P Tla (r, t) 



(3) 



also up to subleading contributions. In this letter 
I only consider color-independent diffusion constants 
Xa — A. Then, for stability, it must be demanded 
J2a "i^aQapn-fj > for all ria > 0. The "temperature" 
variables Tq, measure the difference of the rates of death 
and birth of the species a. Thus the temperatures may be 
positive ore negative. I am interested in the case where all 
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Tq w (up to fluctuation corrections) defining the mul- 
ticritical region. Under these conditions all the species 
live on the border of extinction. 

The equations of motion (|l|) show (absorbing) homoge- 
neous steady state solutions with any combination of ex- 
tinct species. In the infinite volume limit one finds tran- 
sitions to steady states with ria > which are given in 
mean field theory as the solutions of ga/snp = —2Ta 
if tIq > 0. This set of equations in general offers mul- 
tiple different solutions if some or all of the are neg- 
ative. Each of these solutions defines a different homo- 
geneous phase of the active system. In some of these 
phases only one of the species is alive, in the other phases 
several species act cooperatively. In the space spanned 
by the temperature variables Tq, the phases are sep- 
arated by surfaces of first or second order transitions 
where a species becomes extinct. All the critical sur- 
faces meet together in a multicritical point where all 
temperatures are zero. Take e.g. the simplest case of 
a totally symmetric three-species system with g^fs = g 
and Ta < for a, P — 1,2,3. The ternary phase diagram 
is shown in Fig.l. In this case there only exist phases 
where one species is nonzero and the two others are sup- 
pressed. Here stable cooperative phases with more than 
one species are not possible. 

In general, this mean field picture continuously de- 
pends on the coupling constants gap- In the following 
I will show that in spatial dimensions lower than 4, the 
stochastic equations (^^ ) produce universal macroscopic 
steady state equations which correspond to only a dis- 
crete set of choices for the coupling matrix g that I call 
strategies. The stable ones are distinguished by a total 
asymmetry between the members of each pair of species, 
in a mean field picture: gaf3 or gpa is zero. This implies 
that only one of these partners directly influences the 
evolution of the other. These stable strategies with dif- 
ferentiated species permit cooperative phases with more 
than one species nonzero and continuous phase transi- 
tions. All the continuous phase transitions of each strat- 
egy are characterized by the directed percolation critical 
exponents of the onc-spccies Gribov process. 

In the following I will be interested in the behavior of 
correlation- and response functions in the vicinity of the 
multicritical point where all are small. The critical 
properties are extracted using renormalized field theory 
|p^7|--p2| in conjunction with an e-expansion about the up- 
per critical dimension dc = 4. It is then convenient to 
recast the Langevin equation as a dynamical functional 
P3|,p4| 



J [s, s] 



J dtd'^r ^ 



dt + K {to. - 



1^ 1 ^ 



where Sa (r, t) ^ Ua (r, t) are rescaled density variables 
which ensure that gaa = ga- The ia(r,t) denote 
the response fields (the MSR auxiliary variables psf). 
Correlation- and response functions can now be expressed 
as functional averages of monomials of s and s with 
weight exp{— ^}. The responses are defined with re- 
spect to additional local particle sources ha (r, i) > 
in the equation of motion (|l|). By these sources a fur- 
ther term — haSa is added to the integrand of the 
dynamic functional J (^). Coupling constants invariant 
under the rescaling Sa — > CaSa, Sq — > c~^Sa are given 
by /a/3 — gapgi3- The scaling by the usual convenient 
length /i^^ and time scale (A/i^) ^, respectively, shows 

^ Sa ^ M^^^^j /a/3 ~ whcrc E — A — d which sig- 
nals dc = 4 to be the upper critical dimension. A further 
inspection of (^) and a glance at the perturbation ex- 
pansion shows that the functional J includes all relevant 
couplings and is thus renormalizable. 

To determine the renormalizations, one has to cal- 
culate the primitively divergent vertex functions T^^a j 
^aa,a = ^ a,aa i ^ a,ai3 with nonncgativc /i-dimcnsions. 
Here the first and the last part of indices denotes the 
amputated and s-legs, respectively. Dimensional reg- 
ularization and minimal renormalization is used. The 
e-poles are absorbed by multiplicative renormalizations 
of the fields and couplings which I choose in the form 



Aa ^ Aq ^— y \a , Ta ^ — TTTTTT ) 



/a/3 /a/3 = (^) („) rfff "a/3 

/JS Zj-^ Zj^ 



(5) 



(4) 



Here — T {1 + e/2) / (Att)'^^'^ is a convenient constant. 
From the topology of the diagrams of the perturbation 
expansion it is easily seen that the vertex functions do 
not depend on coupling constants with indices other 
than the ones of the vertex functions itself. Therefore 
the renormalization factors Z^"^ with i = (s, r. A) and 
zi"'' := zi"""* that are determined from Ta,a and Taa,a 
are already known from the simple Gribov process. Thus 
the renormalization group flows of the fields and parame- 
ters with exception of the gap with a ^ f3 do not depend 
on the coupling between different colors. It follows at 
once that all the critical exponents are given by the usual 
directed percolation exponents and the fixed point values 
of the "diagonal" renormalized couplings 
all independent from the color a. 

The flow of the couplings between different colors fol- 
lows from the .Z^"'^'* with a ^ (3. The symmetric fixed 



2 



point might be attractive. To disprove that I have calcu- 
lated these renormalizations up to two loop order for the 
special case Aq = A, with the result 

r.(^fi\ 1 ^ 1 r/36 19 



1 3, 4 
e-2l"3 



/ 2 2 



Pa) 



2 A 9 ^ 4 , 

--11 Uaf}Ui3a + l~~4+31n-) Wa/3M 



Collecting together these and the known results for the 
other renormalization factors, one easily gets the flow 
equation for the rcnormalized coupling constants 



^ R 



(6) 



with the Wilson-functions 

/3a/3 = 



1 / N 1 

e + U+ - (Uap + UfJa) - gUa/3W/3o 



97 53 4 

! In - 

128 64 3 
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■a/3 



Uap ■ (7) 



For a — f3, the equation Paa —■ P (ua) = 
yields the nontrivial stable fixed point value = 



{l+( 



169 
288 



53 
144 



In I) e -I- . . .}. For each pair of species 



with colors a ^ f3, fixed point solutions of (j^) with 
at least one zero coupling are either given by Uap* = 
Uf3a* = 0, or by the pair of solutions Uq/j* = (1 + Sap) 
with Eap — —S/3a = ±1 and = -f O i^^)- To ana- 
lyze the case with Uap* and u/ja* both nonzero, one sets 
Uaf3* = u^,+Va[3- To first Order in s one finds a fixed point 
line Vaf3 + vpa — 0. But using the higher order terms of 
(0) , one gets Vap — + O (e^) . Then it can be easily seen 
that Vafi must be zero order for order in the e-expansion. 
This proves that in this case the fixed point is indeed 
symmetric Uq/j* = up^* = *^*- Thus I have found from 
this analysis that for each pair of species independent of 
the other ones three types of fixed points for the couplings 
exist: both couplings vanish (decoupled), or only one of 
the couplings is nonzero (asymmetric), or both couplings 
are nonzero and equal to the diagonal one (symmetric). 

The stability of the various fixed points is analyzed 
by the calculation of the eigenvalues of the matrix of 
derivatives of the Wilson- functions /3a/3, (3 pa- This anal- 
ysis yields as stable fixed points only the asymmetric 
ones. The decoupled fixed points are doubly insta- 
ble, the two instability exponents are equal to ujd — 



(1 



(1 -e 0.0573- 



The sym- 



metric fixed points have one instability exponent lOs = 
(1-3 In I) M 2 . . . = £2 0.0609 + . . .. I especially note 



the very slow crossover from the symmetric to the asym- 
metric fixed point. 

To get a simple qualitative impression of the behavior 
of the dynamic system at its asymmetric fixed points, I 
consider, in the active region with all — Tq, =: gra > 0, 
a rcnormalized mean-field theory of the equations (0J|) 
for spatially homogeneous densities n^and set g^a — g, 
gap = g{l + Sap) and J2a = 



dtUa = \g 



\ ^ Sap flp 

13 



(8) 



The decision which of all the Sap = —£pa are ±1 defines 
a stable strategy. The strategies may be pictured by a 
strategy-diagram. Each corner of such a diagram repre- 
sents a color a. Consider now each pair (a, (3) of corners, 
they are connected by a line with an arrow pointing at a 
if £ap = 1 or vice versa if £ap = — 1. The arrows there- 
fore represent the "pressure" on the reproduction rate 
resulting from one species to another. The diagrams for 
the stable strategies are therefore fully connected by such 
directed lines. For a system with three species e.g. two 
stable strategies are possible, a cyclic and a hierarchic 
one (upper part of Fig. 2). In (H) the influence on the re- 
production rates of the society as a whole is described by 
the common "pressure" p = gn/2. If now Sap — 1 e.g., 
the part of p resulting from the species (3 on the repro- 
duction of a is doubled whereas the part of the common 
pressure resulting from the species a on the reproduc- 
tion of (3 is fully compensated. As a result there exist 
many cooperative phases with more then one species as 
stationary solutions of the mean field equations of mo- 
tion. As an example the phase diagrams for the stable 
strategies of a system of three species is shown in Fig. 2 
(lower part). 

Now the unstable strategies can be pictured as di- 
agrams with disconnected corners for a pair (a, f3) in 
the uncoupled case gap = (then there is no pressure 
from a on (3 and vice versa), and connected corners but 
without direction for a pair (a, P) in the symmetric case 
5a/3 = gpa (then the species a, /? do influence their repro- 
duction rates only by means of the common pressure p). 
It is easily seen that in the fully symmetric case (for three 
species see Fig.l) more then one species can coexist only 
if their temperatures Tq = —gra are all equal. If for one 
color a the parameter is slightly increased, the species 
a suppresses all the others, if is decreased, it becomes 
extinct by the common pressure of all the others. 

Summarily I have shown that also in stochastic mul- 
tispecies models of populations that develops near the 
extinction threshold of all the species and therefore have 
many absorbing states, the critical properties at the mul- 
ticritical point and at all continuous transitions are gov- 
erned by the well known Gribov process (Reggeon field 
theory) exponents (for a previous simulational result on 
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a two-species system which seems to agree with my find- 
ings see |26)). In other regions of the phase diagram of 
course more comphcated critical behavior may arise such 
as multicritical points with different scaling exponents 
(see the recent work of Bassler and Brown |2^]). The 
models considered here have many absorbing states (each 
combination of species may be extincted irrespectively of 
the other ones) and therefore are different from models 
which were considered by Grinstein et al. p8[ where it 
was shown that multispecies systems with one absorbing 
state belong to the Gribov universality class. In addition 
I have shown that the universal properties of interspecies 
correlations and the phase diagram are determined by to- 
tally asymmetric fixed point values of the renormalized 
interspecies coupling constants. This eventually leads to 
a system working cooperatively. It is interesting that 
the asymmetry between the pairs of species seems to be 
the condition for this cooperation near extinction. The 
model considered here is a simple but universal model 
of such a cooperative society and should therefore have 
many applications in all fields of natural science up to 
sociology. 
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FIG. 1. Symmetric strategy of a three-species-system and 
its ternary phase diagram in the plane ri -f T2 -I- ra = — r < 0. 
The numbers a at the corners denote the points where the 
corresponding = —r and the other temperatures are zero. 
Dashed lines indicate discontinuous transitions between the 
one-species phases. 



FIG. 2. Cyclic and hierarchic strategy of a three-species 
system and their ternary phase diagrams. The numbers in- 
side the phase diagrams denote the colors of the living species. 
Sohd lines indicate continuous transitions. 
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